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I. INTRODUCTION 

A renormalizable extension of the Standard Model 
(SM), the vMSM extends the existing SM of particle 
physics by three sterile neutrinos, which are singlets un- 
der the SM gauge group. The i/MSM atteinpts to address 
the various unresolved issues of the SM (for a re- 

view see [5] ) , one of them being dark matter production. 
The lepton asymmetry, produced at the same time as 
the generation of the dark matter sterile neutrinos (at 
temperature of ~ 100 MeV), affects their spectrum and 
number density (see 043)- Upon comparison of theoret- 
ical computation of the abundance of dark matter ster- 
ile neutrino with cosmological and astrophysical observa- 
tions, the lepton asymmetry (AL) is required to be much 
larger than the baryon asymmetry (AB): > 3 x 10^, 
where AB - 10"^°. 

The high value of the lepton-asymmetry can be cre- 
ated provided the masses of the heavier sterile neutri- 
nos in the i^MSM are almost degenerate. Moreover, the 
mass-difference that leads to the requisite value of the 
lepton asymmetry must be much lower than the active 
neutrino mass-difference IT3 | . In an accompanying pa- 
per we discussed the naturalness of the fine-tuning 
required for satisfying the leptogenesis and active neu- 
trino oscillation observations from the perspective of the 
renormalization group (RG) evolution of the i^MSM pa- 
rameters. A complete analysis of the effect of mass- 
difference on leptogenesis, however, requires the compu- 
tation of the physical mass-difference. This is the pur- 
pose of the present work, where we compute the radiative 
corrections to the mass-difference at the one loop level. 
For a review of the generalized on-shell renormalization 
procedure for Majorana neutrino theories, see [isj . 

The paper is organized as follows: in section [TTl we 
review the Lagrangian of the ^MSM, along with the def- 
initions of the different parameters. In section IIIIl we 
describe our formalism of computation of loop correc- 
tions and compute the various contributions. In section 
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^ Different possible ways to modify this requirement may be found 

in [mi. 



ITVl we compute the mass difference with the loop correc- 
tions taken into account, while numerical estimates are 
given in [V] Finally, in section I VII we summarize our 
results. 



II. THE uMSM LAGRANGIAN AND 
RELEVANT MASS-MATRIX 

We use the Lagrangian of the j/MSM in the 
parametrization [3, [l6j . In addition, considering the 
fact that the sterile neutrino Yukawa couplings are much 
smaller than the gauge couplings and neglecting the nu- 
merically much smaller charged lepton Yukawa couplings, 
we choose a basis for the leptonic doublets, in which the 
Lagrangian has the following simple form: 

^uMSM = Cq + A£, (1) 

= CsM + 'N'lid^l^Ni + {hhN2 

1=2,3 

+ f3kN3)4> ~ MN^Ns + h.c, (2) 
A/: = f23hN3^ -^Y. + (3) 

1=2,3 

where Nj are the right handed singlet leptons (/ = 2,3), 
(f) and l2,3 are the Higgs and the lepton doublets respec- 
tively, M is the common mass of the two heavy neutral 
fermions, AM is the diagonal element of the Majorana 
mass matrix, = eijcj)*, M and AM are taken to be real. 
The Yukawa couplings /2, /a can be chosen to be real by 
suitably defining the phases of '3, while the f23 is com- 
plex, with a phase of n. The relation between /2, fs and 
hai introduced in [l,[ll] can be found by comparing eqs. 
© and © with eqs. (2) and (3) of [H. 

We have omitted the dark matter sterile neutrino A^i 
from the Lagrangian as its influence on the problem we 
are interested in is negligibly small f3| . 

We will use unitary gauge ($ = -^^(O v+h)'^) to reduce 
the computational complexity. Here v is the vacuum ex- 
pectation value of the Higgs boson and is taken to be 246 



2 



GeV. For a list of conventions used for the computation, 
see appendix 1X1 

Using the Euler-Lagrange's equation of motion, we get 
the Dirac equation for the system of particles as follows: 
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(4) 



and hence are 8x8 matrices. As usual, mass eigenval- 
ues are given at the one loop level by the poles of the 
propagator. 

Including only one particle irreducible diagrams for 
the computation, it is easy to see that only the fol- 
lowing loop correction matrix elements are relevant: 

T,ss,T,as,^sa and Eaa, where 



{N2RmN2R) {N2RmN3R] 
{N3RmN2R) (iV3fl|E|iV3fl; 



(10) 



and so on. Here the subscripts s, a represent sterile and 
active flavors respectively. 
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At the tree level, the eigenvectors of MIM^ up to first 
order in Yukawa couplings with eigenvalues of M^ up to 
second order in the same are obtained by perturbative 
computation and are listed in the appendix [B] (see fl7|). 

The mass-difference between the sterile neutrino fla- 
vors were shown to be Isl : 



Snitre 



\m 



M ' 



(6) 



B. Loop correction to active neutrino propagators: 
Computation of Eaa 

Considering diagrams that give up to quadratic contri- 
bution in Yukawa coupling, we find that the only possi- 
ble contribution comes from the internal W and Z boson 
loops. The following figure (Fig. [1]) represents the loop- 
correction to the active neutrino propagator. 



w,z 




where 



FIG. 1. Contribution to active neutrino propagator from W,Z 
boson loop. 



hf23v'' + 2MAM. 



(7) 



Performing the one-loop computation using dimen- 
sional regularization, we find that 



The aim of the paper is to compute the first order loop 
corrections to this mass difference. 

Also, the active neutrino mass- difference (Am,y) may 
be solved at the tree- level to be Q: 



{v2l\>^\v2l) = 



(47r) 
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3g^ + / 

24 



(11) 



III. LOOP CORRECTIONS 

A. Propagator 

Consider the propagator for the system of active and 
sterile neutrinos: 



Sf{p) = 



(9) 



where S represents the first order loop correction to the 
mass matrix AA. It is useful to keep in mind that Ai and 
S represent the mass-matrix and the loop corrections of 
the complete system of 2 active and 2 sterile neutrinos. 
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24 
(12) 

Since the active neutrino masses are much smaller than 
the Higgs vacuum expectation value, we can assume that 
<C w^, implying 



1 39' + 9" 
(47r)2 24 



Thus, 



(j/2l|S|i^2l) 



(13) 



(14) 
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Analogous computation shows: 

(i^3l|E|j/3l) = fCi. 



(15) 



In absence of terms contributing to the mixing between 
V2 and 1^3, the lowest order contribution to the loop cor- 
rection matrix may be written as follows: 



San, — 



1 
1 



(16) 



C. Loop Correction to the Sterile Neutrino 
propagators: 
Computation of Ess 

The sterile neutrino propagator receives contributions 
from the internal Higgs loop in addition to W and Z 
boson loops. From here on, we use the intuitive notation 
that T,^,^z,h.t respectively represent the contribution to 
the loop-correction E from W, Z, Higgs-boson loop and 
tadpole graphs. 



3. Contribution from the Higgs boson loop 



FIG. 4. Contribution to sterile neutrino propagator from 
Higgs boson loop. 

The above figure (Fig. |4]) shows the contribution from 
the Higgs loop, which leads to 



{N2R\J:h\N: 



2R) 



■ in ' 
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4-. Matrix Element representing the loop contribution 



Contribution from the W boson loop 



The total loop-correction will be the sum of the afore- 
mentioned contributions. Thus, 




FIG. 2. Contribution to sterile neutrino propagator from W 
boson loop; N, u represent the different possible flavors. 

The contribution coming from the W boson loop is 
represented in Fig. [2] and is given by 



(iV2i?JS»|A^2fl) = 



JL 

{Any 



-1) + 



24p2 



(17) 



where 



Fi = 



-1 



2(47r) 
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(20) 



(21) 
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2. Contribution from Z boson loop 



D. Loop Correction to the Active-Sterile Neutrino 
propagator: Computation of Eas and Esa 




A'b 



The active-sterile neutrino propagator receives similar 
contribution from the W and Z boson loops on the exter- 
nal active neutrino leg and in addition receives contribu- 
tions from tadpoles. We consider only the top quark loop 
among the possible fermion loops in the tadpole graphs. 



FIG. 3. Contribution to sterile neutrino propagator from Z 
boson loop. 



Fig. [3] represents the contribution from the Z-boson 
loop. Performing the computation, we get 



(47r)2 
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4^ Ml ' 



(18) 
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Contribution from W boson 



propagator is represented by Fig. [71 which leads to 

/2 
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FIG. 5. Contribution to active-sterile neutrino propagator 
from W boson loop. 
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4- Matrix Element representing the loop contribution 



The W-boson loop contribution is represented by Fig. 
[5] and it is given by 
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(23) 



Thus, 



^as = (i^l + K2) (^^ = (Xi + i^2)Ct, (26) 



^sa = (i^l + i^2)C. 



(27) 



2. Contribution from Z boson loop 



where 



C = 



/2 
/23 /a 



(28) 
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FIG. 6. Contribution to active-sterile neutrino propagator 
from Z boson loop. 

Fig. [S] represents the contribution from the Z boson 
loop, which leads to 
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3. Contribution from tadpole graphs 



Reasoning as in IIII B[ we arrive at <C . This 
implies 



V 1 3.g2 + £,'2 
"71(4^ 2A 



Nr 




W.Z.h.f 



FIG. 7. Contribution to active-sterile neutrino propagator 
from tadpole graphs. 

The tadpole contribution to the active-sterile neutrino 
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IV. COMPUTATION OF THE 
MASS-DIFFERENCE 

Having obtained the one loop corrections to the prop- 
agator of the system of active and sterile neutrinos, we 
proceed to compute the mass difference between the two 
heavy sterile neutrinos. 

The loop correction can be written in the matrix form 
as follows: 
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In order to obtain the mass-eigenvalues, we find the poles 
of the propagator given in (|9]) i.e. the roots of 



- S = 0, 



(33) 
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Considering only terms up to second order in gauge and 
Yukawa couplings, we arrive at equation of the form: 



'o-y-^-U'^oM^O, ,3., 



r 



1 

V2 



We know that if |e) is an eigenvector of Al},M^, then 

Let 



^ ^ ^ is eigenvector for the matrix ( ^| 



\e) J - yMJ, 

Xi,i — 1 • • • 4 be the eigenvalues for ^ at the tree level. 
Then, applying first order perturbation theory, we obtain 



1 



Thus, 



|W Qi 

[Q2 



^ \i + sx. 



(36) 



(37) 



where, 6Xi represents the change in Xi due to the pertur- 
bation and is a function of . For the on-shell compu- 
tation of the mass-difference, we can safely replace by 
up to second order in Yukawa couplings. 
Denoting the physical mass-difference as SrUpiiys, we 
arrive at a relation: 



^phys 



Srritr 



5m 



loop ; 



(38) 



where ^ j is the unit flavor matrix and Qi, Q2 are 

obtained by multiplying equation p3p with appropriate 
matrices. 

I 
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where neglecting terms of order , 
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Smioop = AM cos aifi + /| + l/asP) + 2/2/23M cos(a + n) {F^ + F2) 
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(40) 



This is central result of this work. As expected, the phys- 
ical mass-difference Sruphys is RG invariant, which can 
be verified using the RG equations for the parameters 
involved. The necessary RG equations are given, for ex- 
ample, in [130 



^ The conventions for the scalar quartic coupling are different in 
the two papers, which must bo kept in mind while performing 
the verification. 



V. NUMERICAL ESTIMATES 

In this section, we provide a numerical estimate of 
the loop contribution to the mass-difference. Evaluating 
equation (j40|) at ^ ^ Mw, we obtain that 



Smioop{Mw) 



2/2/231^' 



(41) 



(47r)2M ■ 

Comparing with our result with ([1]), we see that 

6mioop{Mw) ^ Am^,, (42) 

where Am^ is the active neutrino mass-difference. Thus, 
the loop correction may be absorbed into the Higgs con- 
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densate contribution to the tree-level mass-difference. 
The numerical estimate of the loop correction is taken 
into consideration for the fine-tunings made on the pa- 
rameters of the i/MSM to satisfy the leptogenesis con- 
ditions, as is described in an accompanying paper (see 



VI. CONCLUSIONS 

The generation of low temperature lepton-asymmetry 
has considerable impact on the production of dark mat- 
ter sterile neutrinos. Oscillations or decays of the sin- 
glet neutrinos in the i/MSM can give rise to the requisite 
lepton asymmetry provided their masses are sufficiently 
degenerate. A complete analysis of the problem requires 
computation of the physical mass-difference between the 
heavier neutrinos. In this paper, we computed the loop 
corrections to this mass-difference. On performing the 
computation, we see that the loop correction is of the 
same order as the active neutrino mass-difference. The 
loop-correction is then incorporated into the considera- 
tions for satisfying the leptogenesis conditions, which is 
done in an accompanying paper [l^ . 
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In our convention, the sterile neutrinos are Majorana par- 
ticles with the Dirac spinor written as: 



(A6) 



Thus, PlN = and PrN = Nr. For the active neu- 
trinos, the right handed components do not enter into 
the Lagrangian. So, we can define Dirac spinor for the 
active neutrino as follows: 



(A7) 



Using the properties of Charge conjugation and the com- 
mutation rules for the 75, we get, 

vI^Pru-'.'nI^WPr. (A8) 

Lastly, N'^ = N, in accordance with the fact that N are 
Majorana particles. 



2. Higgs sector 

We use the following form of the Higgs potential: 

= + A((/.V)^ (A9) 

The vacuum expectation value is taken to be w = 246 
GeV. From the above two relations, it can be shown that 
the Higgs mass may be written as: 



(AlO) 



Appendix A: Notations and Conventions 
1. Majorana Fermions 



Appendix B: Mass Eigenstates and Mass 
Eigenvalues of 2 active and 2 sterile neutrino system 

Define the following parametrization: 



^ " j ^ u -1 



a'' = (l,*T),a^ = (l,-^) 



(Al) 



(A2) 



where the Pauli matrices are given by the standard rep- 
resentation 

.^^('l]^.^^('-!]^.^=(l I. (A3) 
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i 



-1 



The left and right projectors are given by: 

^L-^(l + 75),Pfl = ^(l-75). (A4) 



= c^A^, Ct = C-\C^ = -C\ C^-f^C = -7J. (A5) 



./2 = -A; h = -Af; = -eAfe'"; «' = (Bl) 
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This leads to the following form of the mass-matrix: 
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Furthermore, define the parameters: 

fev'^ sin n 



(B3) 



P 



(B4) 



7 



where 
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Then the eigenstates and eigenvalues maybe written as 
follows: 
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